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A generic nonaligned Josephson junction in the presence of an external magnetic field is theoretically con- 
sidered and an unusual flux-dependent current-phase relation (CPR) is revealed. We explain the origin of the 
anomalous CPR via the orbital motion of quasiparticles within a two-dimensional quasiclassical Keldysh-Usadel 
framework. In particular, it is demonstrated that nonaligned Josephson junctions can be utilized to obtain a 
ground-state other than and n, corresponding to a so-called i^-junction, which is tunable via the external mag- 
netic flux. Furthermore, we show that the standard Fraunhofer central peak of the critical supercurrent may be 
inverted into a local minimum solely due to geometrical factors in planar junctions. This yields good consis- 
tency with a recent experimental measurement displaying such type of puzzling feature [R. S. Keizer et al, 
Nature 439, 825 (2006)]. 



PACS numbers: 74.50.+r, 74.45,+c, 74.25.Ha, 74.78.Na 

A sinusoidal current-phase relation (CPR) and Fraunhofer 
response of the critical supercurrent through a s-wave Joseph- 
son contact exposed to an external magnetic field are often 
considered to be standard characteristics of such junctions [1- 
4]. Nevertheless, several theoretical studies have been ded- 
icated to the aim of achieving an experimentally accessible 
situation where the CPR is non-sinusoidal [6-8]. In this case, 
the Josephson ground-state may be characterized by an arbi- 
trary superconducting phase difference ip [6-8], rather than 
the so-called 0- and 7r-states [9]. The first experimental real- 
ization of such a (/^-junction was very recently reported in Ref. 
14. 

Recent studies have also pointed to the fact that the conven- 
tional Fraunhofer pattern in Josephson junctions may be mod- 
ified by the junction geometry [4, 10]. However, there still 
exists experimentally observed magnetic interference profiles 
that remain unsettled in terms of a theoretical explanation of 
the physical origin [11, 12]. In particular, Keizer et al. [11] 
observed an anomalous interference pattern with a local mini- 
mum at zero flux in addition to slowly damped oscillations of 
the critical supercurrent compared to the standard Fraunhofer 
pattern. The setup in Ref. 1 1 consisted of a planar Joseph- 
son junction where superconducting leads were deposited on 
a same side of a half-metallic ferromagnetic strip which was 
fully spin polarized. Fig. 1 A) depicts diagrammatically the 
mentioned experimental setup. To study the system theoret- 
ically, Ref. 13 utilized the Eilenberger formalism in a bal- 
listic planar junction, similar to the setup of Ref. 11, while 
neglecting the orbital motion of quasiparticles. An almost 
^-periodic pattern with non-zero minima of the critical cur- 
rent with respect to external magnetic flux were found due to 
the appearance of second harmonic (sin 2ip, see also Ref. 6). 
However, the inverted interference pattern with a local mini- 
mum at zero flux was not reproduced. 

In this Letter, we consider a generic class of Josephson 
junctions in the presence of an external magnetic field where 
the position of the superconducting leads relative each other 
is not necessarily aligned (see Fig. 1). The obtained results 
are derived without recourse to any ansatz - we have instead 



W1 L W2 L 




FIG. 1: Diagram of considered setups in this Letter. An exter- 
nal magnetic field H (not shown) is applied to the junction in the 
z direction perpendicular to the junction plane. The junction lengths 
and widths are d and W, respectively. A): The planar Josephson 
junction that has experimentally been studied in e.g. Ref. 11. The 
widths of the superconducting leads are assumed to be WIl and 
W— (W1l+W2l)- B): The usual stacked geometry of a Josephson 
junction with displaced superconducting leads. The superconducting 
leads' sizes are WIl and W2r at the top and bottom of junction, re- 
spectively. C): Qualitative view of the quasiparticles' orbital motion 
inside the normal strip subjected to an external magnetic field, which 
is used to describe the origin of the addressed unusual CPR. 

utilized a quasiclassical Keldysh-Usadel technique with the 
numerical approach in Ref. 10 and solved exactly the re- 
sultant linearized equations of motion for the Green's func- 
tion. As our first main result, we unveil that the origin of 
the unexpected interference pattern in the experiment of Ref. 
11 lies within the geometry of the setup. In this way, the 
absence of the standard Fraunhofer pattern, which has not 
been clearly understood, is resolved. In addition to this, we 
demonstrate as our second main result that the CPR in non- 
aligned junctions takes on a very unusual feature: it becomes 
shifted by a term proportional to the external flux <E>, namely 
I(p, <f>) = Io(&) sin((y5 + 0($)) where ip is the superconduct- 
ing phase difference and 9 is a geometry -dependent function. 
Our investigations reveal that the well-known sinusoidal su- 
percurrent and consequently the Fraunhofer pattern manifest 
only in a specific situations. This result is explained in terms 
of the motion of the quasiparticles stemming from the or- 
bital effect. An interesting consequence of the external mag- 
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FIG. 2: Critical supercurrent as a function of external magnetic flux 
<E> through the normal part of junction. The corresponding pair po- 
tential spatial map is given where <p = 0. Throughout the paper we 
have assumed that the junction width is fixed at W=10£s. The first 
and second columns show the critical current I c /Io vs normalized 
external magnetic flux $/$o and corresponding pair potential spa- 
tial maps with various thicknesses d i.e. d=£s and 4£s, respectively. 
Each row indicates different values of WIl and W2l namely, the 
first superconducting lead sizes and the separation of the supercon- 
ducting leads, respectively (see Fig. 1). 



netic flux-shifted superconducting phase-difference is that the 
ground-state of the system may be tuned via the external field 
so that the equilibrium phase difference differs from the con- 
ventional or 7r solutions making a so-called (^-junction. This 
might constitute a simpler alternative to realizing a </?-state 
compared to the array of SFS junctions considered in Ref. 14. 

The starting point for the analysis is the equation of motion 
for the Green's function in the diffusive regime provided by 
the Usadel equation [15]: 

D[d,G(x,y,z)[d,G(x,y,z)]]+i[ep 3 ,G(x,y,z)] = 0, (1) 

where p 3 is 4 x 4 Pauli matrix. Here, e is the particles' energy 
measured from the Fermi level and D is medium diffusive 
constant. In the presence of an external magnetic field H and 
its vector potential A, d = VI — ieA(x, y, z)p 3 provided that 

[d, G(x, y, z)] = VG(i, y, z) - ie[A(x, y, z)p 3 , G(x, y, z)](2) 

[16]. The vector potential is an arbitrary quantity except for 
the restriction V x A = H. We use the Coulomb gauge 
V ■ A = throughout our calculations and assume that the ex- 
ternal magnetic field is oriented in the z direction i.e. H = Hz 
(see Fig. 1). Thus, we may use A = —yHx. In general, the 
Usadel equation should be simultaneously solved along with 
the Maxwell equation V x H = p j in a self-consistent man- 
ner to take into account the influence of screening currents. 
The experimentally relevant scenario is considered where the 
width of the junction W is smaller than the Josephson pen- 
etration length Aj, allowing us to ignore the screening of 
the magnetic field [4, 10, 17]. The Usadel motion equation 
yields a system of nonlinear coupled complex partial differen- 
tial equations that should be supported by suitable boundary 



conditions for studying junctions. In our Josephson system, 
we employ the Kupriyanov-Lukichev boundary conditions at 
N/S interfaces [18] and control the leakage of superconductive 
correlations into the normal strip using an interface parameter 

C; 

C(G(x, y, z)dG(x, y, z)) ■ h = [G B cs(y)> G(x, y, z)], (3) 

in which n is a unit vector denoting the perpendicular direc- 
tion to an interface and tp is the bulk superconducting macro- 
scopic phase. We define ( = Rb/Rf as the ratio between 
the resistance of the barrier region and the resistance in the 
normal sandwiched strip. The bulk solution for the retarded 
Green's function in a s-wave superconductor is given by [16] 
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in which d(e) = arctanh(| A|/£). For a weak proximity effect 
(( ^S> 1), the normal and anomalous Green's functions can be 
approximated by g R ~ 1 and \f R \ <C 1, respectively. The 
current density vector is expressed via the Keldysh block as 

J(K) = J - f deTr{p 3 (G(x,y,z)[d,G(x,y,z)]) K }. (4) 

Here, Jo = N^eD/A and Nq is the density of states at the 
Fermi level. The flux penetrating the junction is given by 
$ = dWH. We also investigate the spatial variation of pair 
potential inside the normal region calculated via: 



U = U Tr{(p x - ip 2 ) / def 3 G K (x,y,z)} 



(5) 



where Uq = —NqA/16 [16]. To study the considered Joseph- 
son junction we use a collocation finite element numerical 
method the same as Ref. 10. The components of approximate 
solution are assumed to be linear combinations of bicubic 
Hermite basis functions satisfying boundary conditions. Ulti- 
mately, the resultant nonsymmetric linear algebraic equations 
are solved via a Jacobi conjugate-gradient method. For more 
details, see Ref. 19. All lengths and energies are normalized 
by the superconducting coherent length £s and superconduct- 
ing gap at absolute zero Ao. The barrier resistance £ is fixed 
at 7 ensuring the validity of weak proximity regime. Temper- 
ature and junction width are T = 0.05T C and W = 10£,s- 
Units also are considered so that ft = ks = 1. 

Figure 2 illustrates the response of the critical Josephson 
current in a planar junction to an external magnetic field as 
shown schematically in Fig. 1 A). Various parameter values 
have been considered in order to make our analysis as gen- 
eral as possible. To do so, we have considered three scenar- 
ios where the superconducting leads have different sizes (first 
row) and where they have equal sizes with a large (second 
row) and small (third row) separation distance. Specifically, 
the third row is relevant with regard to the experiment in Ref. 
1 1 where the size of the electrodes far exceeds the separation 
distance. As seen, in this case the interference pattern ex- 
hibits a local minimum at $ — rather than a maximum as in 
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FIG. 3: Critical supercurrent against external magnetic flux and corresponding pair potential spatial maps of usual Josephson junction with 
displaced superconducting leads including various lead sizes. For the pair potential maps, the superconducting phase difference and external 
magnetic flux are fixed at ip=0 and $=4<i>o> respectively. The junction thickness and width are set to d=2^s and W=10£s, respectively. 




FIG. 4: i) Left column: W1 L = 3£ s , W2 L = 4£ s , WIr = 
3£s, W2 R = 4£ s , ii) middle column: WIl = 6£s, W2 L = 4fs, 
WIr = 0, W2r — 4£s and finally Hi) right column: WIl = 
2£s, W2l = 6£s. The top panels represent the CPRs for various 
values of $/<I>o=0, 3.92, 6.28. The current density spatial maps in 
the bottom row show the results for tp = and <!> = 4$o- The 
superconducting leads' sizes are then set equal at 4£s for all cases as 
schematically are depicted on top of each column. 



the Fraunhofer case, which is fully consistent with the experi- 
mental results in Ref. 1 1 . Whereas it was speculated that this 
minimum might be attributed to a shift in the entire interfer- 
ence curve due to a finite sample magnetization in Ref. 11, it 
is obvious that this is not the case here since the sandwiched 
strip is not ferromagnetic. Moreover, such a shift would make 
the current vs. flux curve manifestly asymmetric (see e.g. Ref. 
20), in contrast to the experimental results of Ref. 1 1 where 
the central minimum is flanked by two large peaks, similar 



to our results. Based on this, it seems reasonable to explain 
the deviation from the standard Fraunhofer pattern as a result 
originating from the combination of a planar geometry with 
the size and separation distance of the superconducting elec- 
trodes. The latter fact is seen by considering the second row of 
Fig. 2 where the separation distance is large compared to the 
superconductors: a Fraunhofer-like pattern emerges, although 
the decay becomes more monotonic as the thickness d of the 
normal strip increases. Even columns in both Figs. 2 and 3 
show the pair potential where the superconducting phase dif- 
ference is zero (p = and external magnetic flux is set to 
$ = 4$o- As seen, the predicted proximity vortices in Refs. 
4 and 10 vanish for the planar junction geometry. However, as 
it will be discussed further below, they reappear in a specific 
condition of stacked geometry (Fig. 1 B)). 

It is instructive to contrast these results with the geometry of 
Fig. 1 B) where the two superconducting leads are connected 
to the normal strip at opposite edges. This is resemblant to 
the experimentally often used stacked geometry. The order of 
frames (critical current and corresponding pair potential spa- 
tial map) are identical to those in Fig. 3 and various lead sizes 
and locations are investigated. It is seen that the location and 
size of both terminals are vital in terms of determining how the 
critical current responds to the external flux. For instance, our 
results reveal that only in specific case where the width of the 
leads' are sufficiently large and connected to opposite edges 
precisely in front of each other does one recover a proximity- 
induced vortex pattern along with the Fraunhofer curve i.e. 

/(</?,$) oc sin $sin^ (6) 

which is a special case corresponding to the scenario of Ref. 
4. The results for the other scenarios in Fig. 3 also show good 
consistency with previous experimental observations [5]. 

Having unveiled the origin of the anomalous inverted 
Fraunhofer response, we now turn to the second main result of 
this paper: the possibility to generate a ^-junction in an SNS 
system with an applied magnetic field. In Fig. 4, we provide 
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the CPR in addition to a spatial map of the current-flow in the 
normal strip for three represented geometries. In i) the leads 
are connected opposite to each other, in ii) they are connected 
antisymmetrically, whereas in Hi) they are connected symmet- 
rically in a planar geometry similar to Ref. 11. It is clear that 
the CPR remains sinusoidal as a function of the superconduct- 
ing phase-difference p in both i) and ii) independent on the 
applied flux. However, case Hi) is qualitatively different. The 
generic form of the CPR is now revealed as: 

=/ ($)simV + e($)) (7) 

in which Iq(&) and 0(<1>) are geometry-dependent functions 
of external magnetic flux as seen in Fig. 4. In fact, Eq. (7) 
holds for all situations considered in Fig. 2 where we have 
demonstrated the CPR is never purely sinusoidal. The stan- 
dard sinusoidal CPR is recovered only for symmetric situa- 
tions relative the induced orbital motion by the external mag- 
netic field, see Fig. 1 C). This observation has a highly in- 
teresting consequence: the anomalous magnetic flux-coupled 
CPR ensures that the ground-state of the system may be tuned 
so that the equilibrium phase difference differs from the con- 
ventional or 7r solutions. Instead, a so-called ^-junction may 
be realized where the ground-state phase difference ip is tun- 
able via the external flux. We therefor arrive at a ground-state 
with Josephson energy ~ 1 — cos(ip + which 

can be controlled by adjusting the applied external magnetic 
field. The idea of a (^-junction via a superconducting phase 
difference shift has been considered previously [7] in the con- 
text of a non-centrosymmetric normal layer with a Rashba 
spin-orbit interaction. However, in our setup the external flux 
is a well-controlled parameter which allows for easy tuning of 
the ground-state, as opposed to controlling a spin-orbit inter- 
action parameter. Moreover, our finding is different from Ref. 
8 where two magnetic junctions, one in 0-state and the other 
in 7r-state with different lengths, are connected in parallel and 
consequently generate an extra cosinusoidal term in addition 
to negative second harmonic. 

What is then the physical origin of this anomalous CPR? 
The answer to this question may be obtained by investigating 
the orbital motion of the qu as ip articles under the influence of 
an external magnetic field inside the normal strip, as seen in 
the current-density maps of Fig. 4 and simple schematic in 
Fig. 1 C). For zero phase difference ip — 0, the external mag- 
netic field induces a current flow where the orbital paths taken 
by the quasiparticles move with the same flux in and out of the 
superconducting regions, in effect no net current flow, only in 
special geometrical configurations. For instance, both in i) 
and ii) the current flow between the superconductors in any 
part of the normal region is seen to have an antisymmetric, 
and thus cancelling, contribution in a different part of the nor- 
mal strip at zero phase difference ip = 0. In contrast, this is no 
longer the case in setup Hi): there is a net flow of current in- 
duced by the orbital response due to the magnetic flux, even at 
(p = 0. To elucidate this clearly in the current-flow, one would 
have to consider the amplitude of the local current as well, but 
the supercurrent-phase curves nevertheless demonstrate that 



this interpretation is correct. In essence, this is a geometry- 
dependent effect since it relies on the positioning of the leads 
relative the induced current-flow via the applied field. Thus, 
it gives rise to the unique possibility to alter the standard CPR 
so that the ground-state of system can be adjusted by tuning 
the external flux. 

To conclude, we have studied the Josephson critical current 
and its response to an external magnetic flux in experimen- 
tally feasible nonaligned junctions. Specifically, a planar ge- 
ometry similar to a recent experiment [11] is considered and 
it is demonstrated that the observed suppression at zero flux 
may stem from the junction geometry rather than any intrin- 
sic magnetization. Moreover, it is shown that a highly un- 
usual supercurrent-phase difference-shift occurs inevitably in 
a class of nonaligned junctions due to an external magnetic 
flux. Its precise form is sensitive to the size and location of 
the superconducting leads. Consequently, this offers a route 
to a tunable junction ground-state. The physical origin of this 
effect is traced back to the orbital motion of the quasiparti- 
cles in the presence of an external field relative the position 
of the superconducting leads. As an interesting consequence, 
this type of Josephson junctions constitute an attainable way 
of realizing the so-called (^-junction experimentally. 
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